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SUMMARY 


This report considers the use of single-degree-of-freedom integrating 
gyros as torque sources for precise control of satellite attitude. Some 
general design criteria are derived and applied to the specific example 
of the Orbiting Astronomical Observatory. The results of the analytical 
design are compared with the results of an analog computer study and also 
with experimental results from a low-friction platform. The steady-state 
and transient behavior of the system, as determined by the analysis, by 
the analog study, and by the experimental platform agreed quite well. 

The results of this study show that systems using integrating gyros 
for precise satellite attitude control can be designed to have a reason- 
ably rapid and well-damped transient response, as well as very small 
steady-state errors. Furthermore, it is shown that the gyros act as 
rate sensors, as well as torque sources, so that no rate stabilization 
networks are required, and when no error sensor is available, the vehicle 
is still rate stabilized. Hence, it is shown that a major advantage of 
a gyro control system is that when the target is occulted, an alternate 
reference is not required . 


INTRODUCTION 


For most satellites some form of attitude control is required, and 
in many cases the control must be fairly precise. This report will 
discuss the use of integrating gyros as torque sources to obtain precise 
satellite attitude control. Similar discussions of systems using motor- 
driven inertia wheels and the earth' s magnetic field are given in 
references 1 and 2. 

The requirements placed on satellite attitude control systems vary 
from one satellite to another; however, it can be stated generally that 
a reference, or tracking line, in the satellite is required to be main- 
tained in alinement with a specified external reference, or line of 
sight, to a specified accuracy. This must be accomplished in the presence 
of disturbing torques and apparent motion of the external reference. 
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Motion of the external reference may "be oscillatory at the orbital 
period, as a result of parallax or velocity aberration, or it may be 
nearly constant, as in the case of an earth-pointing satellite. Dis- 
turbing torques acting on an earth satellite might come from the earth' s 
gravity gradient, the sun's radiation pressure, the earth’s atmospheric 
density, and the earth's magnetic field. 


An earth satellite which is currently being studied by NASA and 
which will serve as an example for gyro system control in this report 
is the Orbiting Astronomical Observatory. The 0A0 will contain one or 
more telescopes, to be used primarily for obtaining photometric and 
spectrographic data from ultraviolet star radiation above the earth' s 
atmosphere. The specific precise control requirements of the 0A0, which A 

will be used as an example in this report, demand that the pointing ^ 

error of the telescope (rigidly attached to the vehicle) be reduced from ^ 

1 minute of arc to less than 0.1 second of arc within 2 or 3 minutes of 3 


time, and the error maintained at less than 0.1 second of arc for one 
orbital period of approximately 100 minutes. Furthermore, roll motion 
about the line of sight must be maintained at less than 1 second of arc 
per second. For solar sighting, the line-of-sLght velocity, due to 
apparent motion of the external reference, wilL have a m ax imum of 0.01 
second of arc per second, predominantly from parallax effects. For 
stellar sighting, the line-of-sight velocity will have a maximum of 
0.005 second of arc per second, predominantly from velocity aberration 
effects. Both of these line-of-sight velocity effects are calculated 
in reference 1. The disturbance torques were assumed to be on the order 
of 100 dyne centimeters. A description of the torque inputs leading to 
this estimate may be found in reference 3- 

In this report a physical explanation of nhe use of a gyro as a 
torque source for a satellite is followed by a derivation and discus- 
sion of the equations for a single-degree-of-f reedom gyro . The over-all 
system equations are then developed, both open loop and closed loop, and 
the system errors are discussed. The equations derived are then normal- 
ized and simplified to allow easy visualization of the system response 
characteristics. A sample system based on the example of the 0A0 is 
designed, and the results of a three-axis analog study are discussed. 
Finally, a comparison is made between theoretical, analog, and 
experimental results for a single-axis system. 


NOTATION 


A 


f 


C 


H 


gyro float angle, radians 

gyro damping constant, dyne cm sec 

gyro angular momentum, dyne cm sec 
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H s stored angular momentum, dyne cm sec 

I-^ moment of Inertia of vehicle, gm cm 2 

If moment of inertia of gyro float, gm cm 2 

k TG over-all gain from sensor to torque generator, (dyne cm) /radian 
s Laplace operator, sec -1 

t time, sec 

t c duration of time during which vehicle is to be controlled, sec 

Uim time during which error signal is limited, sec 

tset settling time, sec 

T c control torque or gyro output torque, dyne cm 

Td disturbance torque exerted on body by external sources, dyne cm 

TlO torque exerted on gyro float by torque generator, dyne cm 

Tjjq uncertainty in the torque exerted on gyro float, dyne cm 

x,y,z axes of a right-hand orthogonal coordinate system 

£ normalized damping ratio 

t time constant, sec 

T n natural period, — , sec 

u)n 

( Pb angle "between tracking line and reference, radians 

cp e error in pointing, radians 

^ e lim effective limit level of error detector, radians 

^LS angle between line of sight and reference, radians 
angular velocity of vehicle, radians/ sec 
angular velocity of gyro float, radians/ sec 
ojls angular velocity of line of sight, radians/sec 
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w n natural frequency of oyer-all control system, radians/ sec 

o vector quantity 

Subscripts 

0 initial value 

01 open loop 

ss steady state 

max maximum allowable value of quantity 

GENERAL SYSTEM DISCUSSION 


A 

k 
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A simplified block diagram of a single-axis control system that 
uses an integrating gyro as a torque source is shown in sketch (a): 


Disturbance 

torque 



Sketch (a) 

The basic command loop consists of a sensor to detect pointing error, 
cp e , an integrating gyro, which produces a control torque, T c , acting on 
the satellite, and finally the satellite itself, which includes a 
rigidly attached telescope. 
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The gyros are also rigidly attached to the satellite, with both the 
gyro float axis and the spin reference axis perpendicular to the axis 
about which it is desired to exert a torque on the satellite. Figure 1 
is a pictorial view of one of the gyros, including the angular momentum 
vector, H, along the spin axis, the float angular velocity vector, wp, 
along the float axis, and the resultant control torque vector, T c , 
along the gyro sensitive axis. The gyro serves both as a source of 
control torque and as a rate stabilizing device. This can be seen as 
follows: If a torque is exerted about the gyro float axis, either by 

the torque generator or by an angular velocity of the vehicle about the 
gyro sensitive axis, the float will move about the gyro float axis at a 
rate, u>f, (assuming a gyro with fluid damping) proportional to the_torque. 
The rate of rotation of the gyro float produces a torque, T c = UfXH, 
about the sensitive axis. When T c is caused by angular velocity of 
the satellite, cjft , then T c will be in a direction to reduce the satellite 
angular velocity resulting in rate stabilization. When T c is caused 
by the torque generator the gyro will again exert a torque on the satellite 
about the sensitive axis. 

The amount of control torque available about the gyro sensitive 
axis will depend on the velocity of the gyro float about the float axis, 
and this in turn, in the steady state, is inversely proportional to the 
damping of the gyro; that is, for a given input torque to the float, if 
the damping is reduced, the steady angular velocity of the float, and 
hence the output torque, will be increased. The system response, however, 
will become more oscillatory; thus, the selection of damping will depend 
on the production of a sufficient amount of control torque without the 
response becoming too oscillatory. 

It should be noted that if it is desired to obtain a constant 
vehicle velocity about the gyro sensitive axis, the gyro float velocity 
must be zero to maintain constant angular momentum (assuming zero disturb- 
ance torque acting on the satellite) . The torque generator must thus 
put out a steady- state torque, T^Qj sufficient to cancel the gyroscopic 
torque, tu^xH, acting on the float. This torque-generator torque will be 
independent of the damping and will react on the vehicle about one of 
the other control axes, producing cross coupling. One possible means 
for eliminating this cross coupling is to use two counterrotating gyros 
for each axis. Then the torque from the torque generator of one gyro 
will cancel that from the other gyro, and there will be no net cross- 
coupling torque on the satellite. However, for the example, where the 
gyros are used to maintain the vehicle essentially inertially fixed in 
space, this cross coupling is found to be negligible. 

There are two disturbances to be considered, as previously mentioned: 
an angular velocity of the line of sight, and an external disturbance 
torque acting on the satellite. Consider first the effect of a constant 
angular velocity of the line of sight. In the steady-state condition, 
the gyro float angle will be constant, as previously noted, and the net 
torque acting on the gyro float must be zero. The torque created by the 
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ang ular velocity of the vehicle about the gyro sensitive axis (to follow 
the line of sight) must be exactly canceled by the torque from the torque 
generator. This latter torque, for the system shown in sketch (a), is 
proportional to the tracking error. 

Next, consider only the effect of external torque on the system. 

To counter this torque, the gyro must have a constant velocity about its 
float axis, which, since the gyro is damped, requires a constant torque 
from the torque generator. Since the satellite is not rotating there 
will be no torque on the float from gyroscopic precession. Again, the 
torque from the torque generator is proportional to the tracking error. 

Thus, to correct for either type of disturbance, there must be a 
steady-state pointing error, and to keep this pointing error small the 
sensor and amplifier gains must be quite large. 


DERIVATION OF SYSTEM EQUATIONS 
Gyro Transfer Function 


The block diagram of an integrating gyro is shown in sketch (b) : 


Amp. and Gyro 

torque gen. Float precession 



Sketch (b) 
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The three summing points shown represent the torque summing action of 
the gyro float, with the net resultant torque being available to accel- 
erate the float, and thus produce an angular velocity of the float, Uf. 
The torque inputs to the gyro about the float axis are from the torque 
generator, Ttq; an uncertainty, Tu G , resulting from gravity unbalance, 
spring torques, etc.j the gyro precession, -uq^H, which is perpendicular 
to and H; and fluid damping, C . Float velocity, ojf, resulting 
from the summation of these torques produces an output control torque, 
T c = w fH, about the sensitive axis. 


The over-all gyro equations obtained from sketch (b) are 


^ ( T ^ + Tu g - Hwb) 
Af 1 + s(l f /c) 

H r T TG + TU G ~ 

T c = HLo f = £ — 

1 C _ i + s(l f /c) J 


( 1 ) 

( 2 ) 


Equation (2) shows that the gyro can have a torque multiplying effect 
(provided H/c > l) , so that the control torque exerted on the body can 
be many times larger than that exerted by the torque generator on the 
gyro float . 


The complete block diagram for single-axis control is shown in 
sketch (c) . 


T U G T d 



Satellite 


Sketch (c) 
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The gyro representation is simplified from that of sketch (b) by use of 
the transfer functions of equations (l) and (2) . Included are the 
external disturbance torque, Ta, acting on the satellite, and the uncer- 
tainty torque, Ty G , acting on the gyro float. The gyro float angle, 

Af, is not used in this application, except as a measure of the amount 
of momentum stored. 

It should be noted that the torque from the torque generator, Ttg> 
the uncertainties, Tu G , and the damping, Ciop, all act on the float and 
react on the case, so that these reactions appear as disturbing torques 
about one of the other satellite axes. However, in the steady state 
their values will be much smaller than the control torque, T c , and so 
will not seriously affect the response of the system. A- 

" 4 

4 

Open-Loop Transfer Functions 3 


The open- loop transfer functions applicable to sketch (e) , assuming 
the loop is broken at point A, are: 

For angular velocity of the vehicle; 


Td|>( s t + f) + (’A +1 W i 


Wb 


Ol 


2 Iblf .IbC 

S, + S h 1 

H 2 H 2 


For gyro float angle; 


^ + Ttc) ?-¥7 


L Ol 


I-blf i b c 

c? *1" S o 

H 2 H 2 


+ 1 


(3) 




Closed-Loop Transfer Fund ions 


If the loop is now closed, the relation 

k TG 

t TG = — ( W LS - u> b) (5) 

can be used to obtain the closed- loop transfer functions. From equations 
( 3 ) and ( 5 ) , the angular velocity equation is found 
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A = 


W LS + ~T 


BK tg 


If 

S ~ + 1 


% 


'G 1 


k TG 


Vf 

HKtg 


4* S' 


P V , JL 

hktg s k tc 


+ l 


( 6 ) 


A 

k 
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Since the chief concern is the magnitude of the error angle, q> g , 
the closed-loop error equation is desired. Using the relationship 

= W LS " sc P e (7) 

and solving for cp e gives the desired equation: 


W LS 


H 

k TG 


I-blf 

H 2 


+ s 


I b C 

H 2 


+ 1 - Ta 


hktg 


If 

s lf + 1 


= 


Iblf 

HKtg 


2 Ib^ H , 

+ s w — + s v — + I 

HKtg k TG 


- Tu G 


k tg 


(8) 


The associated float angle } Af, found by substituting 


sA f H + T d ^ 

Ttg = Ui£ - -■-iJ—) — 


•b 


(9) 


into equation (4) , is 


A„ = 


X b Id ( H \ 1 


sl b 

Tu G HK tg 


• L b x f 

hktg 


I b c 


+ s 


H 


HK TG Ktg 


+ l 


(10) 


Steady-State Considerations 


Consider first the condition when the object star is occulted by 
the earth; that is, Ktg - 0. Applying the final-value theorem for step 
inputs of T[Jq> and Ttg to equation ( 3 ) gives 

^oiss = Id + Iu g 1 + ITG I 


( 11 ) 
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It is Interesting to note that the resultant steady- state satellite 
angular velocity is Independent of the satellite moment of inertia. This 
can be seen physically by realizing that for = 0 the control torque , 

T c , generated by the gyro depends mostly on vehicle angular velocity, 
and when this torque equals the disturbance torque, T d , there will be no 
further change in the vehicle velocity. Furthermore, equation (ll) 
points out a great advantage of the gyro system over other types of 
control systems, namely rate stabilization. Thus, without an error 
detector, the vehicle will assume an angular velocity which varies 
linearly with and is approximately proportional, to the disturbing torque. 
This can be compared to the angular velocity of a vehicle, still with 
no error detector, but also with no gyro aboard. In this latter case 
the external torque produces a constant acceleration, and so the angular 
velocity will increase at a steady rate, that :.s 


UH 


= 2it 


°lss 


( 12 ) 


If equations (ll) and (12) are integrated, the effectiveness of this rate 
stabilization can be determined. After a reasonable period of occultation, 
it is found that for many pratical applications, the error angle with 
the gyro control system is much smaller than tiat for the uncontrolled 
case. In the 0A0 example, the error for the uncontrolled case might be 
much greater than the range of the error detector, while for the gyro- 
controlled case, the error would remain within the range of the error 
detector so that the gyros would eliminate the need for switching to a 
coarse control system following occultation. 


For the occulted condition the steady- state float angle may also 
be found if the inverse transform of equation (4) is simplified for 
large values of t. 


A 

4 

4 

3 




(13) 


This equation shows that the float angle will Lncrease constantly with 
time, and thus that the gyro momentum vector rotates to cancel the effect 
of the external torque. 

For the complete system, with a signal available from the error 
detector, the steady-state equations obtained from equations (8) and (lO) 
are: 


<P e 


C "^G 

e ss KqiQ ^ KipQ. 


H 


■SS 


lb Id 
W LS H “ H t 


(14) 


( 15 ) 
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The sighting error , < Ve S3 ^ reaches a steady-state value which can 

he made small with proper design. The portion of the float angle , Ap , 

s s 

resulting from to^g is a relatively small constant, hut that from the 
disturbance torque increases directly with time. This is to he expected, 
since it shows that the gyro is storing the momentum resulting from the 
external torques. 


NORMALIZED DESIGN 


A 

3 


In the actual design of a control system, the first considerations 
are the amount of disturbance torque expected and the length of the time, 
t c , over which the control system must counter this torque. These two 
factors determine the amount of momentum which the control system must 
he able to store. 


Since torque is the time derivative of momentum, the stored momentum, 
H s , must equal / dt. If a constant disturbance torque is applied for 
the duration of the control time, the maximum stored momentum becomes 
H s = T d t c . In turn, the stored momentum will be the product of gyro 
wheel momentum, H, and the sine of the float angle, Ap. Since the float 
angle must be limited to relatively small angles to avoid excessive 
cross coupling, sine Ap 1 Ap. Hence, 

- ^ Sjnax _ ^d^c 
^^max ^max 

which will allow the determination of H. 



Once a value for H is selected, it seems reasonable to normalize 
all other variables with respect to this value, When this is done, the 
equations previously derived take the following forms: For the open 

loop, from equations ( 3 ) and (4) , 


“b 


T d C ( _ If H 
H C 


H H 1 S 


+ 1 




G 


H 


L TG 

H 


01 *b If lb C 

s ir + s n + 1 


( 17 ) 


A fol 


Tug lb Td 1 

H H / H ~ H s 


s2 ibif + s ibc + 1 
H H H H 


( 18 ) 
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These reduce to the following steady-state equations (from eqs. (ll) 
and (13)): 




01 


ss 


Td C x x Tr JG 

H H H H 


( 19 ) 


Af 

OX SB 



(20) 


For the closed loop, from equations (8) and (lO) , 


H ( 2 lb If . IK , Id C H If H , , V E 

WLS KTG V H H + s H H + y ~ H H K T G \ H C + y ~ H K T G 

s 3 ^b^fJL + s 2 ^ ^ JL - sJL + 1 

H H Kixj H H Ktg K TG 

( 21 ) 


A 

4 

4 

3 


Af = 


J h T d / H A x \3 H 

W LS H " sH VKtg S + 1 ) + H S H Ktq 


S 3 


J b J f H 

IT TT k^g 


+ s 2 


£b C _H_ 
H H K^q. 


+ s 


H 

k TG 


+ 1 


( 22 ) 


And the steady-state -equations (from eqs. (l4) and (15)) become 


9c 


— CO 


'SS 


! ls 


H _^.C _jH_ _ _H_ 

k TG H H K tg H K tg 


Af 


SS 


kJ , 

W LS 1 if 1 


(23) 

(24) 


It should he noted, as mentioned previously, that the term involving 
m ~T c; in equation (24) is usually very small, ard can be neglected in 
comparison with the external torque term. 

Equation (23), which is plotted in figure 2. can be used to determine 
limits on Kg^/ll and (c/H) (h/K^g) by considerirg the values of allowable 
error and external disturbances expected. Since, in general, it is 
desirable to operate without excess forward-locp gain, this figure can 
be used to fix the value of Ktg/H, except that when C/H is selected, 
the product must be checked. The remaining normalized parameters are 
Ib/H, If/H, and C/H. The first of these is known as soon as the vehicle 
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is selected. As will be shown later, it is desirable for the quantity 
If /H 1° l 1 ® small. Its actual size is determined by mechanical consid- 
erations of gyro design. For good system design, c/H should be 
selected on the basis of desired transient performance. This can be 
done more easily if the transfer function is simplified. 

It would be desirable to reduce the form of equation (21) from a 
second-order numerator and third-order denominator to first and second 
orders, respectively. Considering first the denominator, it can be 
shown that if the s 3 term is very small compared to the s 2 term at 
3 — j u n> then, for reasonable values of damping ratio, the s 3 temn may 
be dropped, with negligible effect on the system response. For equation 
(21) this requirement for simplification is stated mathematically as 


£b If _H_ 
H H Ktg 


co. 


n 


< < w n 2 


£b C _H_ 

H H Ktq 


(25) 


which reduces to 


if 

H Wn 


«s 


(26) 


or, alternatively, If/C < < l/co n . Since If/C is the time constant of 
the gyro float response (sketch (c)), it is reasonable to expect that 
If/C will, in fact, be much less than l/co n if the float response time 
is not to degrade the over-all system response. This inequality must, 
of course, be checked when the design is completed. 

The numerator of equation (2l) should now be considered. Since the 
denominator will attenuate greatly all frequencies higher than cj n , it 
seems reasonable to compare the relative sizes of the numerator terras 
at u n . If this is done, the numerator term associated with ojLg, at 
s = J w n> becomes 


I-b I-f 


n 


H H 


TT t TT TT b 1 


lb C 
H H + 


When equation (26) is applied, it can be seen that the first term 
is negligible compared to the second and can be dropped. Further simpli- 
fication of the numerator can be made by comparing the numerator temn 
associated with T^ at s = jw n , 


co. 


If H 


n 


H C + 


1 
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With equation (26) written as 


£f H 
H C 


co. 


n 


< < 1 


it can he seen that the terra associated with T d reduces to unity. Thus 
with the one assumption of equation (26), equation (2l) becomes 


U 1S 


H 

Ktg 


C 

H S 


+ 1 - TT 1 tt 7 ? 


T d C _ JI_ 

h h Ktg 


5* _H_ 

h k tg 


<p* = 


lb C H 

h h Ktg 


+ s 


H 

Ktg 


+ i 


(27) 


This equation is independent of the float inertia. If, and, as would be 
expected from the remarks concerning the inequality of equation (2b) , 
equation ( 27 ) could be obtained directly from sketch (c) by setting 
If = 0. 

Rather common symbols for second-order equations may be used to 
write equation ( 27 ) as 


^e 


W I£ k“- ( Ts + ^ 


T d C H H 

H H% : “ H K tg 


+ 2L s + 
w n 2 w n 


1 


( 28 ) 


vhere 


T 


£b c 

H H 


w 


n 


>_H H k TG 
lb C 


H 


(29) 

(30) 


. _ 1 / H H H 

^ 2 J Ktg c i b 


( 31 ) 


It should be noted that these three quantities are related by 


2£w n r = 1 


(32) 


4=" 4^ > 



CO 4r- 4r- > 
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and further, sub ject to the restrictions of equation (26) , they are 
independent of If/H. 

From dimensionless transient response curves, such as those given 
by Truxal on pages 38 to 4 l (ref. 4 ) , one can see that if 

ToJ n <1 ( 33 ) 

the over-all transient response will have only slightly more overshoot 
than a simple second-order system with no zero. Thus the system can be 
designed for a damping ratio, ^ , in the region of 0.5 to 1.0 (resulting 

in Tco n <1, from equation (32)) with assurance of reasonable transient 
response . 

Selecting a damping ratio of 0.7 (giving Tw n =0.7) and solving 
equation (31) for C/H gives: 

C i H H 

H 2 Kp G I b 13*0 


Changing the desired damping ratio will of course change the numerical 
coefficient in equation ( 34 ) . 

The system design has now been completed, presumably satisfactorily. 

It is desirable to be able to check the behavior of this system when the 
target star is occulted, which is theoretically the same as setting 
K T& = 0- When the inequality of equation (26) is applied to equation (17), 
the following equation is obtained: 


“b 


01 


T d C , . Ttg 

H H + H + H 
s C 
H H 1 


( 35 ) 


This is now effectively a first-order system with a time constant 

x b C 


•oi 


H H “ T 


( 36 ) 


The complete design procedure developed may be summarized as follows: 
(l) select a value for H, using equation (l6) ; (2) select a value for 
k TG/H from figure 2 (a) ; (3) compute the value of C/H required for 

0.7 damping from equation ( 34 ) (or for any other damping ratio from 
equation (31)) ; ( 4 ) check figure 2 (b) to see that the resultant (c/H) 
(H/Ktg) is satisfactory; (5) compute values for r and u n from equations 

(29) and (30); (6) check to see that the inequality of equation (26) is 
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satisfied; (7) check to see that the numerator zero does not deteriorate 
the performance, using equation (33); and (8) examine the open-loop 
performance, using equation (35) • 


ERROR AND TORQUE GENERATOR LIMITING 


In actual practice, the signal from the sensor will "be a voltage 
that will be limited at some relatively small angle, compared to the 
maximum error angle the sensor can detect. A„so, there will be some 
maximum torque, Ttg, available from the torque generator. However, since 
the sensor and torque generator are adjacent , they may be combined as 
shown in sketch (c) . During the time the system is limited, it is 
operating open- loop, and equations (17) through (20) apply. Equation 
(19) shows that if the torque-generator torque (or the error) is limited 
at some maximum value, the vehicle will have a corresponding maximum 
angular velocity. If very large initial errors exist, the gyros will 
accelerate the vehicle to this maximum angular velocity. The vehicle 
will then coast at this velocity until the error becomes less than the 
limiting value. Since the vehicle is thus essentially passive during 
this coast period, the size of the initial error will not affect the 
final portion of the transient after a linear error signal becomes 
available. Thus, system stability is essentially independent of the 
size of the initial error, although the total settling time will be 
lengthened as the initial error is increased. 

An estimate of the time, tq im , during which the error is limited, 
can be made by means of steady-state characte r sties . If the vehicle is 
assumed to rotate at its maximum angular velocity, , until the error 

becomes equal to the limited value, cp e i . , then it can be stated 
mathematically that 


9e “ Te-, • 
B o lim 


"lim 




max 


(37) 


This equation assumes that the time to accelerate to Wb max 
negligible. Equation (19) can be written for this purpose as 

~ t TG _ ^ e lim K T 1 
^max H H 


Thus 


^lim - 


3 e 0 “ ^eqixn 


9 


e lim 


H 

Ktg 


(38) 
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For a second-order system, such as the one represented by equation 
(27), it can be shown that if the damping ratio is about 0.7, the 
settling time from an initial step in error to 2 percent of the initial 
value is comparable to the natural period. If the velocity of the 
vehicle when the error becomes unlimited is neglected and the final error 
is assumed to be about 2 percent of the limiting error, then the linear 
portion of the settling time can be approximated by T n . Thus the total 
settling time may be approximated by 


"^set ~ ^lim + T n (39) 

In this equation T n is too large, since the vehicle velocity at the 
limiting error is assumed to be zero, and ti lm is too small, since the 
initial acceleration is assumed to occur instantaneously. 


SAMPLE DESIGNS 


A set of equations and curves has now been developed which will 
enable a control system to be designed with satisfactory transient 
response and steady-state errors. These curves and equations will now 
be used to design a sample system for the 0A0. The requirements for 
this satellite are as follows: 


<p„ <0.1 second of arc 

e ss — 

t c ~ 105 min 


The satellite moment of inertia is assumed to be 

X-j., = 10 lo dyne cm sec 2 

and the maximum values of disturbances are assumed as 

j oo Tii g = 2.4x10 8 radian/sec ~ 0. 005°/hr 
|t,J = 100 dyne cm 

The necessary gyro constants will be assumed as 


% 


- 1 dyne cm 


If ! -3 

— = 1.4x10 sec 
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To minimize cross coupling, the gyro float angle is limited at 3°. From 
equation (l6) it is found that H = 1.2X10 7 dyie cm sec. Allowing an 
additional margin of safety and using a round number gives 

H = 2X10 r dyne cm sec 


Then 




5X10 8 sec 
5X10 6 sec 


5X10 2 sec 


To meet the steady-state error requirement ((cp es s | < 0.1 second of 

arc) , the following gain and damping requirements are determined from 
figure 2: 


due to ko 


LS 


due to 


due to 


H 


N 

H 


> 0.05 sec 

n — 


>0.1 sec’ 1 

XL — 


HT f ^ 1° sec_1 


For KTf7H = 0.1 sec' 1 , which is the minimum allowable value, 
Ib/H = 5X10 2 sec, and £ = 0.7, it is found from equation (3*0 that 
C/H = 0.01. Checking for the (Ktg/H)(h/c) prc duct gives 


%G H = 0.1 = _ 

H C 0.01 

which is satisfactory. Values for t and 0J n can now be found by using 
equations ( 29 ) and (30) , which give T = 5 sec and w n = 0.14- radian/sec. 
The inequality of equation ( 26 ) becomes 2X10 -4 < < 10“ 2 , which is 
satisfied. 
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The value of the numerator zero is found to he (from eg. (32)) 

Tw n = 0.7. Since this is less than 1, the zero will have little effect 
on the over-all response (as previously discussed for eg. (33)). The 
time constant of the open- loop response is found from eguation (36) as 


A 

4 

4 

3 


^01 — 5 sec 


Thus, the over-all design appears satisfactory, and the table below 
lists the values of the parameters and response characteristics of the 
system. 



KTG 

H 


10 lo dyne cm sec 2 
2xi0 7 dyne cm sec 

l.ifXlO 3 sec 
0.1 sec -1 


i = 0.7 

w n = O.lif radian/sec 
T n = 45 sec 

T oi = 5 sec 


£ = 0.01 

LI 


Egmtions (38) and (39) can be used to estimate the settling time of 
this system. With an initial error of 60 seconds of arc and an assumed 
limit level of 5 seconds of arc , 


t se t = 155 seconds 

This value is not unreasonable, and it can be reduced, if desired, by 
raising the limit level. 


ANALOG SIMULATION 


To check on the analytical work just presented, an analog study was 
made of the gyro control system. The three-dimensional problem was 
simulated; that is, three gyros and three single-axis sensors were 
mounted on the satellite. 
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The gyros were mounted, as shown in sketch (d) , so that the 
sensitive axis of the x gyro was alined along the satellite x axis, 



S = Sensitive axis 
F = Float oxis 
SRA =Spin reference axis 


Sketch (d) 

and the float axis was alined at 45° to the satellite y axis. The 
sensitive axis of the y gyro was along the satellite y axis, and its 
float axis was 45° from the satellite z axis. The z gyro was alined 
in a corresponding fashion. The advantage of this gyro mounting arrange- 
ment (suggested by the Reeves Instrument Corporation) is that the sum of 
the three H vectors is equal to zero. As a result, changes in wheel 
drive frequency will not exert a torque on the satellite. 

Figure 3 shows the analog diagram for the system. Figure 3(a) shows 
one channel of the control system; the other ;wo are identical. Figures 
3(b) and 3(c) show the necessary torque and angular velocity resolutions, 
resulting from the 45° angle between float and satellite axis, and from 
the rotation of the float, Af, about the floa; axis. The broken lines 
in figure 3(a) complete the signal paths when only a single-axis system 
is being considered. 
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Figure 4 shows the response of the system to a step position input 
about one of the axes. It was assumed that all three float angles were 
initially zero and that the error limit was 5 seconds of arc. 1 The 
responses for step inputs about the other two axes had the same charac- 
teristics , with no indication of cross coupling. Runs made with initial 
errors about two or three axes also showed the same response characteris- 
tics, with no cross coupling. 

Figure 5 shows the response to the same input, but this time with 
an initial, very large float angle of the x gyro, A fx , of 15°. The 
input was only in the x axis and, as shown in figure 5> the response 
A ttri. s channel was essentially the same as in figure 4. The cross 

4 coupling produced error and float angular velocity in both the y and z 

4 channels. The responses of these two channels were very similar, although 

3 only the y channel is shown. Additional runs were made in which the 

other gyros had initial conditions, and the responses were very similar 
to those shown in figure 5. 

The effect of torque input was also studied on the analog computer, 
with the results shown in figure 6. Here there was a step input of 
torque of 100 dyne cm, and, as expected, the error angle reached a 
steady-state value of 0.1 second of arc, and the float angle increased 
at a constant rate. The other channels showed no signs of cross coupling. 

Since the cross coupling is small, a study of the effect of varying 
the parameters can be carried out on a single-axis system and the 
predicted response characteristics from the linear analysis should be 
reasonably valid. A sinqplified data analysis of figure 4 indicates a 
system damping ratio, in the linear region, of about 0.7 and a natural 
period of 45 seconds. The settling time, computed from equations (j8) 
and (39) using cp eo — 11 seconds of arc and Teij_ m = 5 seconds of arc, 

is 57 seconds. Thus, the values from the analog simulation agree very well 
with the predicted values. 


EXPERIMENTAL EVALUATION 


Some experimental work was carried out on a low-friction, three- 
degree-of-freedom platform. The platform, shown in figure 7 and described 
in more detail in reference 1, was supported on an air bearing to obtain 
extremely low friction levels. This platform had a considerably smaller 
moment of inertia than the 0A0, and the gyros used on the platform were 
specially modified single-degree-of-freedom gyros. The particular 
constants are 

1 With the value of torque generator gain used, this corresponds 
to a maximum torque-generator torque of 50 dyne cm. 
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It, = 2.4xi0 8 dyne cm sec 2 
H = 9.6x10 s dyne cm se; 
1^ = 1.4xl0 4 dyne cm sec 2 
C = 2.7X10 5 dyne cm sec 


Computing the same normalized design ratios used previously gives 

lb 0 _ 

-g- = 25 sec 

If -3 

-5- = 1.5X10 sec 
il 

% = 0.028 
H 


Experimentally, a gain of K^q/H - 1*7 sec was used to obtain 
reasonable response. From equation (23) it can be shown that the 
following inputs can be tolerated for | cpegg i ^0,1 second of arc: 


^Lsj - 8,4xi0~ 7 radian/sec 
TxJq < 8 dyne cm 
jTd| < 29O dyne cm 


A 

4 

4 

3 


As will be shown later, it was found that the platform had random errors 
of about 1,0 second of arc. From equation (23) , it is found that a torque 
of about 3000 dyne cm would produce a steady- state error of this value. 

The parameters indicated for the system may be used to compute the 
response characteristics. For damping ratio, equation (31) is used to 
obtain 


£ = 0.46 

For natural frequency, equation (30) is used to obtain 


- 1.6 radians/ sec 
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which corresponds to a natural period of 

T n = 3-9 sec 

These values indicate that the response is reasonably damped, but faster 
than necessary. The speed of response cannot be reduced (by reducing 
the gain) however, without causing an increased steady-state error. The 
inequality of equation (26) is satisfied since 

2.4xi0" 3 < < 2.8X10" 2 

The value of T is calculated from equation (29) , giving 

T = 0.7 sec 

and the value of Tw n is found to be 

Tw n =1.1 

This value does not satisfy the inequality of equation (33), so the 
numerator term of equation (28) cannot be completely neglected in 
estimating the response; however, its effect is only to increase the 
initial overshoot by about 15 percent. The time constant, for the open- 
loop system, from equation (36) , is 

Tqi = T = 0.7 sec 


-^fre following table lists the parameters and the system response 
characteristics: 


= 2.4xi0 8 dyne cm sec 2 
H = 9-6xi0 6 dyne cm sec 
If/H = 1.5X10" 3 sec 
C/H = 0.028 
k tc /h =1.7 sec -1 


£ = 0.46 

w n = 1.6 radians/ sec 

T n = 3.9 sec 

T oi = 0.7 sec 


Figure 8(a) shows the experimental transient response of the platform 
to an initial error. Wo limiting occurs in the sensor or torque generator, 
and the natural frequency and damping compare favorably with the computed 
values. It can be seen that accuracy of control is only about 1 second 
of arc. It was felt that this large error was caused by external dis- 
turbances acting on the platform. Figure 8(b) shows the analog simulation 
which used the parameters of the experimental platform and included 
"random" steps of external torque having a magnitude of about 1000 dyne cm. 
It can be seen that the general character of the initial response is 
essentially the same. The amplitude of response of the experimental 
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platform to external disturbances is about three times as large as that 
of the analog system to the analog disturbance. Since the analog dis- 
turbance was 1000 dyne cm, this would indicate an experimental disturb- 
ance of about 3000 dyne cm, which agrees very well with the value of 
disturbance torque previously estimated. The system responses about the 
other axes were similar and no cross coupling between the channels was 
noted. 


CONCLUDING REMARKS 


The results of this investigation have shown that single-degree-of- 
freedom integrating gyros, acting as torque sources, can provide precise 
attitude control of a satellite. The transient behavior of such a system 
is reasonably rapid and well damped. The steady- state errors of the 
system, in response to external disturbances and motion of the line of 
sight, appear to be sufficiently small for many applications. The gyros 
act as rate sensors as well as torque sources , so that no rate stabili- 
zation networks are required, and when no erior sensor is available, the 
vehicle is still rate stabilized. Hence, a major advantage of a gyro 
control system is that when the target is occulted, an alternate reference 
is not required. In addition, because of the torque multiplying effect 
of the gyro, the torque generator and, therefore, the driving amplifier, 
can be smaller than would be the case for inertia wheels, where the 
torque is developed directly by the motor. 


Compared to an inertia wheel system, one disadvantage of the gyro 
system is that only about 5 percent (for a float angle of ± 3 °) of the 
gyro momentum is available for active control; thus a larger or faster 
turning wheel is required. Other disadvantages, which are not mentioned 
in this report, are the additional moving parts and weight required for 
the gyro float and case, which may have to be temperature controlled, 
requiring additional power. These advantages and disadvantages, along 
with total system weight and power, must be considered when a system for 
an actual satellite is being designed. 


A 

4 

4 

3 
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(a.) oj t . q and Tuq as disturbances. 

Figure 2.- Gain requirements vs. steady-state error 
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(b) T d as disturbance. 


Figure 2.- Concluded. 



Figure 3.- Analog diagram of three-axis control system 
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(b) Torque resolution. 


Figure 3*- Continued. 
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(c) Angular velocity resolution. 
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Figure 3*- Concluded. 
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Figure 7.- Satellite attitude control test platform with gyros. 
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(a) Experimental system response. 
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(b) Analog system response. 


Figure 8.- Response of platform to an initial step position input 
accompanied by disturbance torques. 
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